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ABSTRACT 


The  Caucasus-Caspian  region  is  part  of  the  Alpine-Himalayan  collision  belt  and  is  an  area  of  complex  structure 
accompanied  by  large  variations  in  seismic  wave  velocities.  Using  data  from  29  new  broadband  seismic  stations  in 
the  region  as  well  as  data  from  a  temporary  (1999-2001)  deployment  in  eastern  Turkey,  a  unified  velocity  structure 
is  developed  using  teleseismic  receiver  functions  and  surface  waves.  Joint  inversion  of  surface  wave  group 
dispersion  curves  generated  from  ambient  noise  with  receiver  functions  show  that  crustal  thickness  varies  from  34  to 
52  km  in  the  region.  The  thickest  crust  is  in  Lesser  Caucasus  and  the  thinnest  is  in  the  Arabian  Plate.  Thin  crust  is 
also  observed  near  the  Caspian.  The  lithospheric  mantle  in  the  Greater  Caucasus  and  the  Kura  depression  is  faster 
than  the  Anatolian  Plateau  and  Lesser  Caucasus.  This  possibly  indicates  the  presence  of  cold  lithosphere.  The  lower 
crust  is  slowest  in  the  northeastern  part  of  the  Anatolian  Plateau  where  Holocene  volcanoes  are  located. 

Fundamental  mode  Rayleigh  wave  phase  velocities  are  determined  at  periods  between  20  and  145  seconds.  We 
observe  a  relatively  high-velocity  zone  located  in  the  upper  mantle  under  the  Kura  basin  and  the  western  part  of 
Caspian  Sea  that  is  continuous  to  the  Moho.  The  images  show  very  low  velocities  beneath  the  eastern  Anatolian 
plateau  imply  ing  the  existence  of  a  partially  molten  asthenospheric  material  underlying  a  very  thin  lithosphere. 

Using  a  two-station  method,  both  Lg  and  Pg  attenuation  is  measured  and  tomograpically  inverted  to  yield 
attenuation  maps.  Efficient  Lg  propagation  is  observed  throughout  much  of  the  Arabian  plate.  Moderate  Lg  Q  is 
observed  in  the  Lesser  Caucasus  and  Kura  Basin  while  low  Lg  Q  is  observed  in  the  East  Anatolian  plateau.  Pg 
shows  highly  variable  propagation  throughout  the  region. 
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OBJECTIVES 

The  objective  of  the  research  is  to  develop  3D  seismic  velocity  structure  and  regional  phase  attenuation  models  for 
the  Caucasus-Caspian  region  (Figure  1).  Specifically,  the  purpose  of  this  work  was  to: 

•  Create  upper  mantle/crustal  velocity  models  using  data  from  new  stations  in  the  region. 

•  Construct  detailed  maps  of  regional  phase  attenuation. 

•  Compare  and  validate  results  and  models  using  the  various  algorithms  as  well  as  independent  datasets 


The  work  consisted  of  four  main  tasks:  data  collection,  regional  phase  analysis,  crustal  and  upper  mantle  velocity 
determination,  and  model  validation. 


RESEARCH  ACCOMPLISHED 

Data  eolleetion.  Waveform  data  was  collected  from  33  stations  (Figure  1)  in  the  region  as  part  of  collaborative 
effort  with  the  individual  seismic  surveys.  The  primary  focus  was  collecting  continuous  data  from  2006-2008, 
although  waveform  data  from  other  years  was  collected  as  available.  All  data  has  been  provided  to  the  Knowledge 
Base  at  Lawrence  Livermore  National  Laboratory. 


1UUy  G*Orgto 

Figure  1.  (left)  Area  of  study  showing  broadband  stations.  Broad-band  stations  include  STS-2,  ESP3T, 
ESP3ESD  and  CMG40T  type  instruments.(right)  Example  receiver  functions  stacked  by  baek- 
azimuth  for  selected  Azerbaijan  stations.  Global  map  shows  recorded  events  (red  dots)  and 
location  of  stations  (yellow  star).  Small  red  circle  shows  position  of  Saatly  8  km  deep  w  ell  used 
for  model  validation. 


Surfaee  w  ave  group  velocity  and  receiver  function  joint  inversion.  To  develop  a  comprehensive  model  of  crustal 
and  upper  mantle  velocities,  joint  inversion  of  receiver  functions  and  surface  waves  was  applied.  Crustal  receiver 
functions  are  sensitive  to  velocity  discontinuities  while  surface  wave  dispersion  is  controlled  by  average  S-wave 
velocity  structure  (e.g.  Julia  et  al.,  2000).  Therefore,  combining  the  two  methods  yields  a  robust  estimate  of  crustal 
and  upper  mantle  velocities.  Receiver  functions  using  events  at  distances  between  30  and  90  degrees  were  calculated 
for  all  stations  using  iterative  deconvolution  at  Gaussian  widths  of  1.0,  1 .5,  and  2.5  seconds.  Each  receiver  function 
(Figure  1 )  was  examined  for  consistency  and  signals  with  excessive  amplitudes  on  the  tangential  component  or 
grossly  different  from  the  stacked  average  for  each  station  was  discarded.  Coverage  was  best  for  back-azimuths  to 
the  east  and  sparse  to  the  west.  Initially,  the  receiver  functions  were  modeled  independently  to  investigate  quality. 
The  slant  stacking  of  Zhu  and  Kanamori  (2000)  was  also  used  to  estimate  Moho  depths.  In  general,  the  stations  near 
the  Caspian  show  poor  quality  receiver  functions  with  indications  of  considerable  crustal  scattering  and 
reverberations  due  to  the  thick  sediments.  These  were  difficult  to  model.  Some  indications  of  3D  structure  were 
observed  for  stations  in  the  Greater  Caucasus  as  the  receiver  functions  showed  clear  variations  with  back  azimuth. 
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correlations  with  respect  to  GNI  (filtered  10-  50  s) 
red  -  paths  to  east 


■HO  -JO*  100  0  K»  »  M  JOQ 


lag  (seconds) 

Figure  2.  (top)  Ray  paths  used  in  ambient  noise  surface  wave  correlation,  (bottom,  left)  Correlations 
of  stations  with  respect  to  station  GNI.  Stations  marked  in  red  show  ray  paths  from  west  and 
stations  in  yellow  are  east  of  GNI.  (bottom,  right)  Dispersion  curves  for  stations  pairs  located 
in  the  Anatolia  plateau  and  across  the  South  Caspian. 


Surface  wave  group  velocity  dispersion  measurements  were  made  using  both  event-based  and  ambient  noise 
correlation  methods.  Love  and  Rayleigh  wave  group  velocity  dispersion  curves  were  estimated  for  over  1,500 
waveforms  at  periods  of  7-100  see.  Ambient  noise  correlation  was  applied  to  all  possible  station  pairs  (-400)  to 
obtain  Rayleigh  wave  group  velocities.  All  dispersion  curves  were  picked  manually  and  then  included  in  the 
global/rcgional  tomographic  inversion  of  Pasyanos  et  al.  (2005).  We  then  extracted  the  dispersion  curves  from  the 
tomography  maps  of  surfaec  waves  and  combined  them  with  stacked  receiver  functions. 


Horizontal  depth  slices  of  1 1 , 35  and  87  km  are  shown  in  Figure  3  as  well  as  an  estimated  Moho  thickness  map.  The 
thick  sediments  of  Kura  Basin  are  evident  in  upper  crust  with  low  velocities  averaging  Vs=2.8  km/s.  The  eastern 
part  of  the  Greater  Caucasus  shows  similar  low  velocities  in  the  upper  crust.  The  slowest  lower  crustal  velocities  arc 
observed  in  the  northeastern  Anatolian  plateau  and  Lesser  Caucasus  region  where  considerable 
Neogene/Quaternary,  Holocene  volcanic  activity  has  occurred.  Faster  lower  crustal  velocities  occur  at  the  edge  of 
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the  Caspian  and  in  the  Arabian  plate.  The  upper  mantle  appears  distinctly  slow  under  the  Anatolian  plateau  and  to 
the  southeast.  Indications  of  faster  upper  mantle  velocities  occur  under  the  Greater  Caucasus.  It  was  noted  that 
inversion  for  Love  and  Rayleigh  waves  provided  slightly  different  results,  possibly  due  to  anisotropy. 


Figure  3.  A,  B)  Crustal  S  velocities  estimated  from  joint  inversion  at  depths  of  1 1  (left)  and  35.5  km  (right).  C) 
Upper  mantle  S  velocities  at  a  depth  of  87  km.  D)  Moho  depth  inferred  from  velocity  profiles  based 
from  joint  inversion.  Black  lines  mark  international  borders  and  coastlines;  red  lines  show  major 
faults  and  suture  zones.  Black  dots  are  Holocene  volcanics. 


Surface  wave  phase  velocity  inversion.  To  investigate  variations  in  upper  mantle  velocities,  surface  wave  phase 
velocities  measured  from  earthquakes  were  also  inverted.  Data  from  an  earlier  temporary  deployment  of  29  stations 
was  used  as  well  as  the  newer  data.  Teleseismic  earthquakes  with  a  magnitude  greater  than  5.8  were  selected  and  the 
waveforms  were  evaluated  for  high  signal  to  noise  at  the  longer  periods  (50  and  125  seconds).  Unfortunately,  many 
of  the  Azerbaijan  stations  showed  high  noise  levels  at  the  longer  periods  even  on  the  vertical  component,  possibly 
due  to  barometric  variations.  This  restricted  the  data  set  to  20  events.  The  data  for  each  event  were  filtered  to  create 
13  frequency  bands  with  corresponding  centered  periods  ranging  from  20  to  145  seconds  and  dispersion  curves  were 
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estimated  from  these  results.  Initially,  the  Rayleigh  wave  phase  dispersion  was  inverted  to  estimate  ID  velocities. 
The  results  were  topographically  inverted  for  2D  structure  on  a  50  km  grid  using  two  methods:  Forsyth  et  al.  (1998) 
and  Yang  and  Forsyth  (2006).  The  two  methods  are  similar  but  the  Yang  and  Forsyth  technique  attempts  to 
compensate  for  scattering  caused  by  local  heterogeneities  and  therefore  should  improve  spatial  resolution  (Figure  4). 
The  results  at  shorter  wavelengths  were  more  robust  than  for  the  longer  wavelengths  especially  for  the  area  covered 
by  the  Azeri  stations,  possibly  due  to  excessive  noise  at  longer  wavelengths.  Higher  upper  mantle  velocities  were 
observed  under  Anatolia  and  lower  velocities  under  the  Caucasus/Caspian. 


0  00  0  02  0.01  0  04  0  05  0C6  0  07  0  08  1  00 

Figure  4.  (top)  Preliminary  results  from  S  wave  velocities  inferred  from  Rayleigh  wave  phase 

dispersion  at  periods  of  50  (left)  and  125  seconds  (right).  The  method  of  Yang  and  Forsyth 
(2006)  was  used.  The  lower  plots  show  the  estimated  errors  in  units  of  velocity  with 
triangles  denoting  station  locations.  Due  to  high  noise  at  longer  periods,  results  from  the 
Caucasus  at  longer  periods  are  not  well  constrained. 
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Regional  phase  analysis.  The  purpose  of  the  regional  phase  analysis  was  to  extend  regional  phase  attenuation  maps 
(e.g.,  Al-Lazki  et  al.,  2003;  Gok  ct  al.,  2003)  throughout  Anatolia  and  the  Caucasus  (Figure  5).  The  mapping  was 
performed  by  first  measuring  attenuation  between  stations  and  then  inverting  to  solve  for  the  average  Q  in  discrete 
spatial  cells.  The  direct  two-station  method  was  used  (e.g.,  Zor  et  al.,  2007;  Xie  et  al.,  2004)  to  estimate  attenuation 
for  each  path.  This  method  relies  on  spectral  ratios  between  pairs  of  station  aligned  with  the  ray  paths  from  each 
event.  Phases  were  identified  manually  to  avoid  problems  caused  by  variations  in  crustal  velocities.  Spectra  for  each 
phase  were  estimated  after  windowing.  Station  pairs  for  each  event  were  selected  using  a  criteria  based  on  intcr- 
station  distance  and  alignment  with  event  (±15°).  The  spectral  ratio  between  two  stations  was  calculated  for 
appropriate  station  pairs  after  correction  for  distance  (Xie  et  al.,  2004)  and  geometrical  spreading.  Q  and  q 
(frequency  dependence)  were  then  estimated  from  the  spectral  ratios  using  linear  regression.  This  Q  and  q  represents 
the  attenuation  along  a  line  between  the  two  stations.  These  attenuation  measurements  were  added  to  an  existing 
Middle  East  dataset  which  was  inverted  to  create  a  regional  map  (Figure  5).  The  dataset  includes  approximately 
3000  station  pairs.  Separate  inversions  were  conducted  for  both  Lg  and  Pg.  Checkerboard  tests  were  performed  to 
evaluate  resolution  and  both  ART  and  LSQR  inversion  algorithms  were  tested  to  evaluate  the  effect  of  possible 
numerical  artifacts.  Results  were  similar  for  both  algorithms.  Raypath  coverage  is  good  in  eastern  Anatolia  but  is 
sparser  towards  the  Caucasus  and  Caspian. 
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The  1  Hz  Lg  results  show  a  broad  zone  of  high  attenuation  (low  Q)  extending  roughly  east-west  north  of  the  Arabian 
plate  from  Anatolia  to  the  Caspian  (Figure  5).  Considerable  spatial  variation  exists  which  likely  reflects  the 
complicated  tectonics.  Lg  appears  to  propagate  well  in  the  Arabian  plate  but  is  dramatically  attenuated  in  the  Lesser 
Caucasus.  This  may  be  due  to  the  volcanism  and  the  possibility  of  partial  melt  in  the  crust.  The  Kura  basin  and  the 
extreme  eastern  Caucasus  show  moderate  to  low  attenuation,  which  may  reflect  the  influence  of  the  Eurasia  plate.  It 
is  clear  from  these  results  that  the  Kura  basin  is  underlain  by  continental  crust  and  is  distinct  from  the  South 
Caspian.  The  variations  in  Lg  attenuation  appear  to  result  from  a  combination  of  both  crustal  heterogeneity  and 
changes  in  crustal  thickness.  While  Moho  depth  increases  by  5  to  10  km  at  the  northern  edge  of  the  Arabian  Plate 
the  variation  in  depth  from  Anatolia  to  the  Lesser  Caucasus  is  fairly  small  and  yet  relatively  low  Lg  Q  is  observed. 
Volcanic  regions  are  often  observed  to  possess  low  Lg  Q  possibly  due  to  high  attenuation  of  shear  waves  by  partial 
melt  in  the  crust.  Pg  reveals  a  distinctly  different  pattern.  The  transition  between  the  Arabian  plate  and  Anatolian 
plateau  is  defined  by  an  abrupt  change  in  Pg  Q  east  of  about  38  E  longitude  but  a  clear  zone  of  high  Q  exists  in  the 
Lesser  Caucasus.  The  zone  of  low  Q  along  the  Caspian  and  Kura  basin  may  be  caused  by  the  abrupt  thickening  in 
sediments  and  dipping  Moho,  which  decreases  the  efficiency  of  the  crustal  P  waveguide. 


0  1'JU  2UC  3UU  4U'J 

Lg  Q  at  1  Hz 

Figure  5.  Attenuation  maps  for  Lg  (left)  and  Pg  (right). 


Model  validation.  The  model  was  tested  in  several  ways:  comparison  with  alternate  datasets  such  as  reflection 
seismic  and  gravity,  earthquake  hypocenters  comparison,  and  with  waveform  modeling.  Seismic  reflection 
refraction  data  shows  clearly  the  deep  sedimentary  layers  in  the  South  Caspian  and  coastline  (Knapp  et  al.,  2004; 
Mangino  and  Priestley,  1998).  Gravity  data  also  suggests  deep  basement,  which  closely  matches  the  low  velocities 
for  the  area  observed  by  the  joint  inversion.  The  Moho  depths  inferred  by  refraction  data  are  also  similar.  2D 
modeling  of  gravity  data  through  the  center  of  the  Kura  basin  shows  a  clear  gravity  high  that  runs  roughly  north- 
south  and  parallel  to  the  Caspian  coast  and  delineates  the  Kura  from  the  South  Caspian.  The  existence  of  this 
basement  high  is  confirmed  by  the  Saatly  deep  well,  which  penetrated  to  a  depth  of  9  km  through  a  succession  of 
sediment  and  then  volcanoclastic  layers.  To  first  order,  alternate  datasets  are  consistent  with  our  results,  both  for 
crustal  and  mantle  velocities  and  with  phase  propagation. 
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A  second  set  of  tests  was  conducted  with 
earthquake  locations.  This  effort  was 
focused  in  the  Azerbaijan  region.  The  3D 
model  was  averaged  to  produce  an 
average  Id  model  Hypoeenter  locations 
were  calculated  using  this  model  and 
compared  with  global  PDE  (Preliminary 
Determination  of  Epicenters)  locations. 

Significant  differences  in  location  were 
observed  (up  to  30  km)  (Figure  6).  A 
major  factor  is  the  exceptionally  thick 
sedimentary  package  in  Kura  Basin  and 
nearby  South  Caspian  Basin.  Waveforms 
for  events  with  known  foeal  mechanisms 
were  also  calculated  using  this  model. 

CONCLUSIONS  AND 

RECOMMENDATIONS 

An  improved  erustal  and  upper  mantle 
velocity  and  attenuation  model  has  been 
developed  for  the  Caucasus/Caspian 
region  and  Anatolian  plateau.  The 
Caspian  and  Kura  basins  show 
pronounced  low  velocities  in  the  upper 
crust  but  slightly  elevated  lower  crustal 
velocities.  Crustal  thickness  varies  from 
about  50  km  in  the  Lesser  Caucasus  to 
35  km  in  the  Arabian  plate.  The  upper 
mantles  display  strong  heterogeneities  in 
the  region.  It  is  slow  under  the  Anatolian  plateau  and  slightly  fast  under  the  Greater  Caucasus.  These  are  likely 
related  to  slab  remnants  and  ongoing  volcanism. 

Striking  variations  in  regional  phase  propagation  and  velocities  are  mapped.  Lg  is  highly  attenuated  not  only  in  the 
South  Caspian  and  Black  Sea  basins  but  also  in  the  Anatolian  plateau,  possibly  due  to  a  combination  of  crustal 
properties  under  the  volcanic  areas  and  varying  crustal  thickness.  A  clear  change  in  Lg  attenuation  occurs  at  the 
northern  edge  of  the  Arabian  plate.  A  narrow  band  of  moderate  Lg  attenuation  exists  in  the  Kura  basin  and  along  the 
Caspian.  Pg  also  shows  variations  in  attenuation. 

Due  to  the  strong  spatial  variations  in  both  erustal  and  mantle  properties,  regional  seismic  waves  should  be  expected 
to  vary  significantly  both  in  travel-time  and  amplitude  in  this  region.  More  seismic  stations  are  scheduled  to  be 
installed  in  the  region  by  the  various  surveys  and  this  data  would  be  useful  in  increasing  the  resolution  of  these 
results. 
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Figure  6.  Comparison  of  global  earthquake  locations  (stars, 
PDE)  with  locations  determined  using  the  local 
network  (circles)  with  the  new  velocity  model.  Red  is 
shallower  than  25  km,  green  between  25  and  50,  and 
blue  is  deeper  than  50  km.  Lines  connect  different 
locations  of  the  same  event. 
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